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Chimera states occur in networks of coupled oscillators, and are characterized by having some 
fraction of the oscillators perfectly synchronized, while the remainder are desynchronized. Most 
chimera states have been observed in networks of phase oscillators with coupling via a sinusoidal 
function of phase differences, and it is only for such networks that any analysis has been performed. 
Here we present the first analysis of chimera states in a network of planar oscillators, each of which 
is described by both an amplitude and a phase. We find that as the attractivity of the underlying 
periodic orbit is reduced chimeras are destroyed in saddle-node bifurcations, and supercritical Ifopf 
and homoclinic bifurcations of chimeras also occur. 
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Networks of coupled oscillators and their synchroniza- 
tion properties have been studied for many years [l|, Q . 
One particular class of interest involves phase oscillators, 
where each oscillator is described by a single angular 
variable 0, Hj. The use of such phase models is justi- 
fied when the attraction to an underlying limit cycle is 
"strong" relative to the effects of other oscillators in the 
network [H, i, i]. Recently a number of investigators 
have studied "chimera" states in networks of phase os- 
cillators in which some fraction of the oscillators 
synchronize while the remainder run freely, even though 
the oscillators may be indentical. Early analyses of these 
states d, 0, [13, IT2l - [l7| used a self-consistency argument 
which can be traced back to Kuramotoja to show ex- 

[Ti[r 



istence of chimeras. Later work d, | 



used the 
SI to de- 



remarkable ansatz of Ott and Antonsen 
rive differential equations governing the evolution of or- 
der parameters of the systems under study, allowing one 
to determine the stability of chimera states and the bi- 
furcations they may undergo. 

It has long been known that networks of identical phase 
oscillators, coupled through a sinusoidal function of phase 
differences, have non-generic behaviour jT8l. l22l424| . Most 
chimera states have been observed in such idealized net- 
works, and in order to determine whether chimeras might 
be observed in real physical systems one should investi- 
gate their robustness with respect to, for example, het- 
erogeneity in intrinsic frequencies, or variations in oscil- 
lator arnplitude. The first issue has already been ad- 
dressed [aQ, and here we investigate the second. 

Several authors have observed chimeras in networks of 
oscillators described by more than one variable [1, 0, [H, 
HgI . [25I [26j . so they are known to exist, but these au- 
thors have either provided no analysis, or have reduced 
their (identical) oscillators to phase oscillators in order 
to analyse their dynamics using the approaches men- 
tioned above. In this Letter we give the first analysis of a 
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chimera state in a network of planar oscillators in which 
the reduction to phase oscillators is not performed. 
The model we consider is 
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for j = -I- 1, . . . 2N, where Xj E C, and oj, e, a, ji and v 
are real parameters. 

These equations describe a pair of populations of N 
Stuart-Landau oscillators with all-to-all coupling within 
each population of strength /x, and all-to-all coupling be- 
tween the two populations of strength v. Such oscillators 
are related to the normal form of a Hopf bifurcation, and 
are a specific example of A — a; oscillators [1, ^, 2J} ■ Such 
a pair of coupled populations of oscillators has been stud- 
ied by several authors [III, HI, HI] , and can be thought 
of as the simplest "network of networks" that one could 
study. 

Defining Xj — Vjc"^^^ , Eq. ([1]) can be written 
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and Eq. ([2]) can be written as a similar pair of equa- 
tions. From Eq. ^ we see that as e 0, the rate of 
attraction to the hmit cycle Vj — 1 becomes infinite, 
and Eq. (U) reduces to Equation (1) of jjjj (after a re- 
definition of w), i.e. our system reduces to a previously- 
studied network of phase oscillators. We will investigate 
the dynamics of Ql-© when e 7^ 0. By allowing the 
radius r to vary, we expect a wider variety of behaviour 
than that seen in networks of phase oscillators; for ex- 
ample, oscillator death and chaos [s^l. For comparison 
with previous results we define (3 = 7r/2 — a and we let 
/i = {l+A)/2, V — {l~A)/2, where A is a parameter 

Firstly, we show a chimera state for (IT])-(l2]); see Fig. [1] 
Panel (a) shows a snapshot of all 9j at an arbitrary time. 
We see that population two (with N -\- \ < j < 2N) has 
completely synchronised (all Vj w 1.0019), while oscilla- 
tors in population one (with 1 < j < N) remain incoher- 
ent. Panel (b) shows that oscillators in population one 
lie on a closed curve (a slight distortion of the unit cir- 
cle) in the complex plane. Panel (c) shows the angular 
density of the oscillators in population one. It is non- 
uniform, i.e. these oscillators are not completely incoher- 
ent, and it was the dynamics of this density that Abrams 
et al. [ll| studied, using the parametrisation of Ott and 
Antonsen In this chimera state the oscillators in 

population two have a constant angular velocity and ra- 
dius, and the distributions in panels (b) and (c) of Fig. [T] 
remain stationary. It is worth noting that the chimera 
state shown in Fig. [I] is attracting, i.e. nearby states are 
attracted to it, unlike the corresponding chimera states 
in networks of identical phase oscillators which are neu- 
trally stable [III, [l^. Allowing both the radius and 
the phase of the oscillators to vary seems to eliminate 
the non-generic behaviour seen in networks of identical, 
sinusoidally-coupled phase oscillators, in the same way 
that making the oscillators non-identical does [1, Q . 

We briefiy digress to analyse the chimera state shown 



in Fig. [T] in the limit e 



0, 



IVj. Let 9 J 



e 



for + 1 < j < 2N and move to a coordinate frame 
rotating with angular velocity fl in which Q is constant. 
Using rotational invariance, set = 0. Then, from the 
equation for population two. 
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and (using Eq. ([5])) each oscillator in population 1 satis- 
fies 
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where S 



= A^"^Ef=iSin(6'fc -a) and C = 
N~^^'^^-^cos{ek ~ a). In the limit N 00, S and 
C are constant and can be replaced by the expected val- 
ues of sin {e — a) and cos {e — a) respectively, calculated 
using the angular density, p{e), which is proportional to 
the reciprocal of the velocity, de/dt 0,[l3|- Thus chimera 
states are described by the simultaneous solution of 
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FIG. 1: (color online) A chimera state for HI)-®, (a): A 
snapshot of the 9j. (b): Tj as a function of 9j (relative to 
^jv+i) for j — 1, . . . N. (c): The density of the 6j's, relative 
to 9n+i, for j = 1, ... A''. Parameters: A*" — 500, lj — 0,e = 
0.05, 13 = 0.08, A = 0.2, 5 = -0.1. 



and 
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and K is a. normalization fac- 
= 1. Following solutions of 



where p{e) = K(de/dt)-^ 
tor such that J^^ p{e)de 
Eqs. (IZl)-® as parameters are varied one can find re- 
gions of parameter space in which chimera states exist, 
in agreement with the results of Abrams et al. [ll| (re- 
sults not shown). Eq. ^ can be interpreted as saying 
that in a chimera state, each oscillator in population one 
follows a periodic orbit, and is nonlinear ly driven by its 
own mean field. This effect is known to be capable of 
destroying completely synchronous behaviour [3l|. We 
now analyse the chimera state in (HJ-lHl) for e 7^ using 
a similar argument, showing that it can be described by 
a single complex number. 

Let Xj =Y ioi N + 1 < j < 2N and go to a rotating 
coordinate frame such that Y is constant in this frame. 
Rotate the frame so that Y is real and positive. Then 
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from Eq. ([2]) we have 

= iico -n)Y + €-^ {1- {l + Sei)Y'^}Y 

+ e-"' ( ^lY + i^x) (9) 



where X = N ^ X^feLi -^k, and each oscillator in popu- 
lation one satisfies 
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Given X, the real part of Eq. (|9]) can be solved for Y , 
and the imaginary part of Eq. (j9l) can be used to show 
that each oscillator in population one satisfies 



dX 
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+ {l-(l + 5ez)|X|2} X + e-'^'l^iX + vY] (11) 



i.e. each oscillator in population one is driven in a non- 
linear way by the mean field of population one. Thus our 
self-consistency equation, i.e. the analogue of ©-([H]), is 
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where X(t; X) is a periodic solution of Eq. ([TT|) with pe- 
riod T{X). The main difference between Eqs. (HI)-® and 
Eq. is that X{t] X) must be found by numerically 
integrating Eq. (fTTj) to find a periodic solution, whereas 
the periodic solution of Eq. ^ need not be found — 
only the density, p(0), proportional to the reciprocal of 
the angular velocity, is needed. 

Having found a solution of Eq. (jl2p , it can be numeri- 
cally continued as parameters are varied. Typical results 
are shown in Fig. [2] where we vary e. We see that for 
these parameter values the solution of Eq. ([T^ can be 
continued to e « 0.109, where it appears to undergo a 
saddle-node bifurcation. For e small, points on the lower 
branch in panels (a) -(e) correspond to the stable chimera 
known to exist |11| when e = 0, while the upper branch 
corresponds to the saddle chimera. A typical solution of 
Eq. ((ni) is shown in Fig. H (f). 

The saddle-node bifurcation seen in Fig. [5] can be fol- 
lowed as a second parameter, say (5, is varied. The result 
is shown in Fig. |3] (dashed curve). We see that as 5 is 
increased, the range of values of e for which a chimera 
state exists also increases. However, the curve of saddle- 
node bifurcations in Fig. |3] relates only to the existence 
of chimeras (found through a self-consistency argument 
similar to that of Kuramoto 5] ) not their stability. Nu- 
merical simulations of Eqs. (HJ-© show that a stable sta- 
tionary chimera which exists to the right of the dashed 
curve in Fig. |3]can undergo a supercritical Hopf bifurca- 
tion as parameters are varied, leading to a "breathing" 
chimera [ll|. These oscillatory states then seem 
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FIG. 2; (color online) The solution of Eq. gH). (a): Re(X); 
(b): Im(X); (c): Y; (d): Q. and (e): T{X), as functions of e. 
(f): Real and imaginary parts of the self-consistent solution 
of Eq. (jlip for parameter values shown with a circle in panels 
(a)-(e). Parameters: 13 = 0.08, A = 0.2, 5 = -0.01. 



to be destroyed in a homoclinic bifurcation as parame- 
ters are further varied. Numerically determined curves 
of Hopf and homoclinic bifurcations are shown in Fig. |31 
These curves are conjectured to terminate at a Takens- 
Bogdanov bifurcation on the curve of saddle-node bifur- 
cations, which seems to be the generic arrangement for 
chimera states [1, [ll|, [l^ . Varying A or /3 rather than 
5 results in a similar arrangement of saddle-node, Hopf 
and homoclinic bifurcation curves (results not shown). 

To the left of the dashed curve in Fig. |3] and above 
the curve of homoclinic bifurcations, the perfectly syn- 
chronous state {Xj — Xk Vj, k) is stable. Despite the 
radii of our oscillators being able to vary, we have not 
been able to find oscillator death or more exotic dynam- 
ics by varying parameters. Perhaps this is not too sur- 
prising, since non-identical oscillators (which we have not 
considered here) and strong coupling relative to the at- 
traction to the limit cycle (i.e. the opposite limit from 
that considered here) seem to be required to observe os- 
cillator death [30, S3 ■ 
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FIG. 3: (color online) Bifurcation curves in the 5 — e plane for 
chimera solutions of ([l|-(l2|. Hopf and homoclinic bifurcations 
were found by direct simulation of ([I])-®. A — 0.2, /? — 
0.08, N = 500. 

In principle, the stability of the chimera states studied 
here, and thus the location of the Hopf bifurcation seen in 
Fig. [3l could be determined using the ideas presented in 
Sec. 6 of Matthews et al. [30]. However, a difficulty arises 
because we do not have an analytic expression for the 
chimera state around which to linearise — the density, 
p(r, 9), of oscillators in population one can only be found 
indirectly by numerically solving Eq. . (Note that the 
stability or otherwise of the periodic solution of Eq. ([TT]) 
that we find is not related to the stability of the chimera 



state. Solving Eq. ([TT|) is just a convenient way of finding 
the invariant density for population one.) 

For chimeras to be observable in a physical system they 
must be generic, and not only occur in networks of iden- 
tical phase oscillators with all-to-all coupling via a sinu- 
soidal function of phase differences, which are known to 
have unusual properties [3, H^-Hjl- Their persistence 
when phase oscillators are made non-identical has been 
charcterised previously [8,, ^9,] , and in this Letter we have 
shown that chimeras also persist (within limits) when 
both the amplitude and phase of the oscillators are al- 
lowed to vary. 

One caveat is that the system studied here has all-to-all 
coupling, both within and between populations. It would 
be interesting to determine whether this is necessary in 
order to observe chimeras. Indeed, this raises a more 
general question as to which network topologies support 
chimeras. Also, the system ([D)-(ll]) is invariant under the 
phase shift Xj e'^'^Xj Vj, where 7 is a real constant. 
This seems to be the reason that, in a chimera state, the 
synchronous population undergoes uniform rotation at 
fixed radius in the complex plane, and we can describe the 
incoherent population as having a stationary distribution 
in a uniformly rotating coordinate frame. It would be of 
interest to study chimeras in networks for which this is 
not the case. Addressing these two issues would help 
determine the general robustness of chimeras, and thus 
the likelihood of them having relevance to the physical 
world. 

I thank Steve Strogatz for correspondence which in- 
spired the work presented here. 
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